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Abstract 



We study and simulate M = 2 supersymmetric Wess-Zumino models in one and two di- 
mensions. For any choice of the lattice derivative, the theories can be made manifestly 
supersymmetric by adding appropriate improvement terms corresponding to discretiza- 
tions of surface integrals. In one dimension, our simulations show that a model with 
the Wilson derivative and the Stratonovitch prescription for this discretization leads 
to far better results at finite lattice spacing than other models with Wilson fermions 
considered in the literature. In particular, we check that fermionic and bosonic masses 
coincide and the unbroken Ward identities are fulfilled to high accuracy. Equally good 
results for the effective masses can be obtained in a model with the SLAC derivative 
(even without improvement terms). 

In two dimensions we introduce a non-standard Wilson term in such a way that the 
discretization errors of the kinetic terms are only of order 0(a 2 ). Masses extracted 
from the corresponding manifestly supersymmetric model prove to approach their con- 
tinuum values much quicker than those from a model containing the standard Wilson 
term. Again, a comparable enhancement can be achieved in a theory using the SLAC 
derivative. 
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1 Introduction 



Supersymmetry is an important ingredient of modern high energy physics beyond the stan- 
dard model; since boson masses are protected by supersymmetry in such theories with chiral 
fermions, it helps to reduce the hierarchy and fine-tuning problems drastically, and within 
grand unified theories, it leads to predictions of the proton life-time in agreement with 
present day experimental bounds. As low energy physics is manifestly not supersymmet- 
ric, this symmetry has to be broken at some energy scale. However, non-renormalization 
theorems in four dimensions ensure that tree level supersymmetric theories preserve super- 
symmetry at any finite order of perturbation theory; therefore, supersymmetry has to be 
broken non-perturbatively [I]. 

The lattice formulation of quantum field theories provides a systematic tool to investigate 
non-perturbative problems. In the case of supersymmetric field theories, their formulation is 
hampered by the fact that the supersymmetry algebra closes on the generator of infinitesimal 
translations (2j. Since Poincare symmetry is explicitly broken by the discretization, one 
is tempted to modify the supersymmetry algebra so as to close on discrete translations. 
However, as lattice derivatives do not satisfy the Leibniz rule, supersymmetric actions for 
interacting theories will in general not be invariant under such lattice supersymmetries. The 
violation of the Leibniz rule is an 0(a) effect, and supersymmetry naively will therefore be 
restored in the continuum limit. In the case of Poincare symmetry, the discrete remnants of 
the symmetry on the lattice are sufficient to prohibit the appearance of relevant operators in 
the effective action which are invariant only under a subset of the Poincare group and require 
fine tuning of their coefficients in order to arrive at an invariant continuum limit. In generic 
lattice formulations, there are no discrete remnants of supersymmetry transformations on 
the lattice; in such theories, supersymmetry in the continuum limit can only be achieved by 
appropriately fine-tuning the bare couplings of all supersymmetry-breaking counterterms |3j. 
An additional complication in the formulation of supersymmetric theories on the lattice is 
the fermion doubling problem. Local and translationally invariant hermitean Dirac opera- 
tors on the lattice automatically describe fermions of both chiralities [U [5]; the fermionic 
extra degrees of freedom are usually not paired with bosonic modes and so lead to supersym- 
metry breaking. Generic prescriptions eliminating these extra fermionic modes also break 
supersymmetry. 

As a simple supersymmetric theory, the Wess-Zumino model in two dimensions has been 
the subject of intensive analytic and numerical investigations. Early attempts included 
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the choice of the nonlocal SLAC derivative (thereby avoiding the doubling problem) in the 
Af = 1 version of this model [6l [7j [8j. This was motivated by the idea that the Fourier 
transform of the lattice theory should coincide with that of the continuum theory. A Hamil- 
tonian approach where the SLAC derivative minimizes supersymmetry-breaking artifacts 
introduced by non-antisymmetric lattice derivatives was discussed in [9]. Alternatively, a 
local action with Wilson fermions was constructed at the cost of a nonlocal supersymme- 
try variation [l0]if| simulations of this model (T7J indicate that this theory with a cubic 
superpotential indeed features non-perturbative supersymmetry breaking. 
In order to manifestly preserve some subalgebra of the Af = 1 supersymmetry algebra on the 
lattice, the above discussion suggests to choose a subalgebra independent of the momentum 
operator. Thus, supercharges Q + and Q- can be defined using (non)local derivatives on a 
spatial lattice. The choice of a continuous time then allows for a Hamiltonian defined by 
H := Q 2 + = Q 2 _ which commutes with the supercharges [HI [19] and automatically contains 
a Wilson term. This strategy can be generalized to the Af = 2 model on a spatial lattice; 
an analysis shows that a subalgebra admitting an 0(2) R-symmetry can be preserved |20j. 
Simulations for the Af = 1 (see [2T]) and the A" = 2 model [22J have been done using the 
local Hamiltonian Monte-Carlo methods. Further simulations based on the Green function 
Monte-Carlo method for A" = 1 indicate that supersymmetry is unbroken for a quartic 
superpotential and has a broken and an unbroken phase for cubic superpotentials [23l [24] . 
With a discrete time-coordinate, one has to resort to action- rather than Hamiltonian-based 
approaches. The perfect action approach was pursued for the free A" = 1 theory in [25J; 
a generalization for interacting theories seems however problematic. A treatment of the 
A" = 2 Wess-Zumino model in the Dirac-Kahler formalism preserves a scalar supersymmetry 
on the lattice but leads to non-conjugate transformations of the complex scalar field and 
its conjugate. This enlarges the space of states and presumably renders the theory non- 
unitary [10J. A related [26j idea avoiding these problems goes back to the idea of Nicolai [27] 
that (in this case) scalar supersymmetric field theories admit new bosonic variables with a 
Jacobian cancelling the determinant from integrating out the fermions, in terms of which the 

3 A similar approach with staggered fermions leads to problems in the continuum limit [H]. - In the 
four-dimensional model with TV = 1 supersymmetry with Ginsparg- Wilson fermions, lattice chiral symmetry 
is incompatible with Yukawa couplings |12l 113): however, the theory can be regularized by supersymmetric 
higher derivative corrections, which leads to a supersymmetric continuum limit within perturbation the- 
ory [TJ]. It could be shown that supersymmetric Ward identities are satisfied up to order 0(g 2 ) in the 
coupling constant |15| I16j. 
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bosonic action is purely Gaussian. In general, the new variables are complicated nonlocal 
functions of the original ones; for J\f = 2 supersymmetry, however, local Nicolai variables 
can be expected [28]. Starting from a discretized form of the Nicolai variables found in the 
continuum, one simply defines the bosonic part of the action as the sum over squares of the 
Nicolai variables; the fermionic part of the action is adjusted so that the determinant of the 
fermion matrix still cancels the Jacobian. This leads to an action manifestly preserving part 
of the continuum supersymmetry on the lattice, and which contains improvement terms in 
addition to the naive latticization of the continuum action [28], [20]. As stated in [29, 30j, the 
breaking of some of the continuum supersymmetries can in this framework be traced back to 
the fact that the improvement terms are not compatible with reflection positivity. However, 
the violation of Osterwalder-Schrader positivity is an 0(a 2 g p h ys ) effect and thus should be 
negligible at least at weak coupling [31]. An analysis of the perturbation series shows that 
these terms in an off-shell formulation of the theory with a cubic superpotential and Wilson 
fermions lead to tadpole diagrams which diverge linearly in the continuum limit [32]. A 
cancellation between these would-be central charge terms and the naive discretization of 
the continuum Hamiltonian has been suggested as a solution for the conundrum raised 
in [29] that the lattice result of the number of zero-modes of the Dirac operator seems to 
differ vastly from the continuum answer [9]. 

The possibility to introduce Nicolai variables is closely related to the fact that the (2, 2) 
Wess-Zumino model is a topological theory of Witten type, i.e., the action is of the form 
S = QA for a scalar supercharge Q, once auxiliary fields are introduced [33], [34]. This 
formulation manifestly preserves Q-supersymmetry on the lattice and therefore guarantees 
that the theory remains supersymmetric in the limit of vanishing lattice spacing. Lattice 
theories of this type in various dimensions and with different degrees of supersymmetry 
have been classified in [35]. It turns out that in this formulation some remnant of the U{l)y 
R-symmetry which is left unbroken by the superpotential at the classical level but broken 
by the Wilson terms is restored in the continuum limit even non-perturbatively [HI]; this 
also indicates that at least for the cubic superpotential under consideration supersymmetry 
is not broken nonperturbatively. 

At a fixed point in space, the Wess-Zumino model reduces to supersymmetric quantum 
mechanics. A naive discretization of the action with a Wilson term for just the fermion 
leads to differing masses for fermions and bosons in the continuum limit [36]. This can be 
traced back [37J to 1-loop contributions (of UV degree 0) of the fermion doublers to the boson 
propagator which have to be cancelled by an appropriate counterterm that was neglected 
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in [36]. As soon as the action is amended by a Wilson term also for the boson [36], an 
additional diagram involving the boson doublers cancels the finite correction of the fermion 
doublers, and no counterterms are required to achieve the desired continuum limit. In 
this case, fermion and boson masses agree in the limit of vanishing lattice spacing, and 
supersymmetric Ward identities are fulfilled to great accuracy. 

As expected for a theory with J\f = 2 supersymmetry, supersymmetric quantum mechanics 
can be formulated in terms of local Nicolai variables [28j. However, this action differs from 
the amended form of [36] by a further interaction term which becomes an integral over a 
total derivative in the continuum limit. 

In this paper, we scrutinize six different lattice actions for supersymmetric quantum mechan- 
ics based on Wilson and SLAC fermions with and without such improvement terms. The 
three models without improvement terms are not manifestly supersymmetric in the pres- 
ence of interactions and differ by bosonic terms which become irrelevant in the continuum 
limit as well as by the choice of the lattice derivative; the other three models with manifest 
supersymmetry differ in the prescription for the evaluation of the improvement term and 
also in the choice of the lattice derivative — they can be constructed from three different 
Nicolai variables. We compare the effective masses for interacting theories with a quartic 
superpotential and analyze the Ward identities of the broken and unbroken supersymmetries 
at various lattice sizes and couplings. The central observation here is that the manifestly 
supersymmetric theory with the so-called Stratonovich prescription [38] for the evaluation 
of the improvement term (which is the discretization of a continuum surface integral) leads 
to far better results than the model with an Ito prescription at finite lattice spacing. The 
mass extraction for the models with SLAC derivative is at first sight hampered by an os- 
cillating behavior for nearby insertions in the bosonic and fermionic two-point functions; 
however, this Gibbs phenomenon is under good analytical control and can be softened by 
the application of an appropriate filter. With the help of this "optimal filter", the results 
even surpass those of the model with Stratonovich prescription. 

Along the way, we show for which superpotentials and derivatives one can guarantee posi- 
tivity of the fermion determinants. In those cases where the determinant can be computed 
exactly, we analyze under which circumstances they converge to the correct continuum re- 
sults; these exact results are crucial for our simulations. Again, the model with Stratonovich 
prescription is ahead of the one with Ito prescription; the former leads to a determinant with 
the correct continuum limit, whereas the latter differs from it by a factor which depends on 
the superpotential. 
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Furthermore, we study three different manifestly supersymmetric discretizations of the Af=2 
Wess-Zumino model in two dimensions. Instead of trying to generalize the Stratonovich pre- 
scription to two dimensions, we introduce a non-standard Wilson term (corresponding to 
an imaginary Wilson parameter in the holomorphic superpotential) in such a way that the 
discretization errors for the eigenvalues of the free (bosonic and fermionic) kinetic operators 
are only of order 0(a 2 ) instead of order 0(a) for the standard Wilson term. In the simula- 
tions, we study the effect of the resulting violation of reflection positivity and compare the 
results with those of the model with SLAC fermions. Due to calculational constraints, we 
have to restrict the computations to smaller and intermediate values of the coupling. 
As a theoretical background, we show that the discretized Wess-Zumino model in two di- 
mensions with the SLAC derivative has a renormalizable continuum limit. 
The paper is organized as follows: In section [2], we introduce the quantum mechanical mod- 
els on the lattice, with and without improvement terms, and discuss which behavior of the 
interacting theories at finite lattice spacings can be gathered from their (non-)invariance 
under supersymmetry transformations in the free case. In subsection 12.41 we give details 
about the positivity of the fermion determinants and derive their respective continuum lim- 
its. The results of the effective masses and the Ward identities can be found in section [31 In 
section HI we discuss the discretizations of the Af = 2 Wess-Zumino model in two dimensions 
and present the results of the mass extractions. Section [5] covers the algorithmic aspects of 
our simulations including a derivation of the Gibbs phenomenon for the SLAC correlators 
which justifies the application of the filter for the mass extraction in the quantum mechani- 
cal model. Finally, section [6] contains the proof that the Wess-Zumino model on the lattice 
with the SLAC derivative is renormalizable to first order in perturbation theory; a technical 
part of this proof is completed in appendix [A] 

2 Supersymmetric quantum mechanics 

In this section we introduce the action for supersymmetric quantum mechanics in a lan- 
guage which can be easily generalized later on to the two-dimensional Euclidean Wess- 
Zumino model. In subsections 12.2.11 and 12.3.11 we present six different lattice versions of the 
continuum theory with Euclidean action 

S con t= fdr^ + ^ + ^ + ^W"^, where W'{<j>) = (1) 
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The continuum model is invariant under two supersymmetry transformations, 

5M<f> = eij), 5^ij = -e(j> + W), 5^ip = 0, 

with anticommuting parameters £ and £. The lattice approximations considered below differ 
by the choice of the lattice derivative and/or the discretization prescription for a continuum 
surface term; in subsection 12.41 we argue why three of them lead to far better approximations 
to the continuum theory. The results of our simulations for a quartic superpotential W(<j)) = 
y0 2 + |0 4 (with positive m, g) are discussed in section [3l 



2.1 Lattice models 

We start from a one-dimensional periodic time lattice A with real bosonic variables 4> x and 
two sets of real Gratemann variables i/; x , %[) x on the lattice sites x G A = {1, . . . ,N}. The 
integral and continuum derivative in (pp) are replaced by a Riemann sum a ^2 an d a lattice 
derivative d, where a denotes the lattice constant. Two different antisymmetric lattice 
derivatives will be used in what follows. These are the ultralocal derivative 

9xy = 2^(^+l»» ~ fix-l,y) (3) 



with doublen 4 



and the nonlocal SLAC derivative without doublers, which for an odd number 



of lattice sites takes the form [39], [9] 

d± c v = - '- -—t 1 ttttt and 8^ = 0. (4) 

a sm(n{x -y)/N) 



If we allow for non-antisymmetric derivatives then we may add a multiple of the symmetric 
lattice Laplacian 

A x y = —^{Sx+l^y — 25 xy + 8 X -\ t y) (5) 

to d to get rid of the doublers. In this way, we obtain a one-parameter family of ultralocal 
derivatives with Wilson term, 

8-^arA, -1 < r < 1, (6) 



4 It should be noted that the symmetric combination of forward and backward derivatives leading to ([3]) 
yields an antisymmetric matrix (d xy ). 
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interpolating between the forward (or right) derivative for r = — 1 and the backward (or 
left) derivative for r = 1. As will be argued in subsection 12.41 for the quartic superpotential 
W(4>) = tt0 2 + |0 4 with positive parameters m and g to be considered below we shall need 
the latter with matrix elements 

®xy ~ &xy ~ Ij^-xy — ~{$x,y ~ $x— l,y)- (7) 

The backward derivative is free of doublers and not antisymmetric. For periodic boundary 
conditions the derivative operators <9 slac , d and d h are all given by circulant matrices which 
commute with each other. 



2.2 Lattice models without improvement 

A straightforward discretization of the continuum action (pQ) would be 

Snaive = ~ J] + W%) + a £ (3* + W xy ) (8) 



2 

x,y 



ie lattice derivatives d as well as for 
so that the theory is free of fermion 



where we scrutinize below three possibilities for t 
the terms W x , W xy derived from the superpotentia' 
doublers. None of these models is supersymmetric under the discretization of any of the 
continuum supersymmetries (j2J), 



6W<f> x = eip x , 5^4> x = -E((d<P) x + W x ), 5«^ = 0, 

5® <\> x = $ x e, 6®# x = 0, 8^^ x = {{d<j)) x - W x )e, 1 1 

however, at least for a free theory both supersymmetries are realized in two of the models 
with antisymmetric matrices (d xy ). Thus, we might expect a better approximation to the 
continuum theory for these models. In fact, it will turn out in section [3] that this behavior 
pertains to the interacting case, e. g., the masses extracted from these two models are much 
closer to their continuum values than those from the third theory. Truly supersymmetric 
(improved) lattice models will be considered in section I2~3l 

2.2.1 The unimproved models in detail 

(i) Naive lattice model with Wilson fermions 

The most naive discretization is given by the action (JSj) with an additional Wilson term 



5 In general, W x is not equal to W'(4> x ). The definition for each model can be found below. 
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shifting the derivative d as explained in (JJJ), i.e., 

Sf} = S^e, 9 = <9 b , = W (</»,) and W xy = W"((f) x )5 xy . (10) 

The Wilson term removes fermionic as well as bosonic doublers. This action has no su- 
persymmetries at all, and bosonic and fermionic excitations have different masses in the 
continuum limit in the presence of interactions. Even the free model has no exact super- 
symmetry; this can be traced back to the fact that the derivative <9 b is not antisymmetric. 

(%%) Naive lattice model with shifted superpotential 

An alternative way to remove the fermion doublers employed by Golterman and Petcher 
and later Catterall and Gregory [10], [36] is to use the (unshifted) antisymmetric matrix (d xy ) 
and add a Wilson term to the superpotential, 

sgW naive , d = d, W x = ~(A<f>) x + W'((f> x ), W xy = -^A xy + W"(<j )x )6 xy . (11) 

It should be noted that (as compared to (fTUjl ) only the bosonic terms are changed. This 
model is only supersymmetric without interaction, i. e, for W'((j) x ) = m(j) x . In the inter- 
acting case all susy Ward identities are violated. The breaking is equally strong for both 
supersymmetries. 

(Hi) Naive lattice model with SLAC derivative 

Naively, one might expect the supersymmetry breaking effect in the naive lattice action ((HJ) 
with SLAC derivative (111) to be of the same magnitude as with backward derivative d h . 
Surprisingly enough, this is not so; it will turn out that the mass extraction from the model 

SS = Snaive, d = d^ C , W x = W'(<j> x \ W xy = W"(<P x )S xy (12) 

is about as good as for the improved actions considered below. This is again related to the 
fact that the derivative is antisymmetric such that the free model with SLAC derivative 
admits both supersymmetries (in contrast to the model with backward derivative). 

2.3 Lattice models with improvement 

In order to preserve one of the two lattice supersymmetries in ((9j) for interacting theories 
the naive discretization (fSj) should be amended by extra terms which turn into surface terms 
in the continuum limit and reinstall part of the continuum supersymmetry on the lattice. 
Such invariant lattice models may be constructed with the help of a Nicolai map i— > £(0) 
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of the bosonic variables. In terms of the Nicolai variables £ x (4>) the improved lattice actions 
take the simple form 

Ssusy = ^^ix{4>) 2 + aJ2^x^-^y (13) 

x x,y y 

This is a discretization of the most general supersymmetric action in terms of a real bosonic 
variable and two real Graftmann variables on the circle [27]. It is easily seen to be invariant 
under the first type of transformations 

5 {l) <px = S {l) i>x = = 0. (14) 

For the particular choice £ x (4>) = [d<p) x + W x the supersymmetric action f TT3l becomes 

5 S usy = < ^^{{d<P) x + W x f + aJ2^ (d X y + W xy )^ y 

x x,y 

= S naive + aJ2WM(d<P) x (15) 

X 

and the supersymmetry transformation ffT4l is identical to 5^ in Qj. So the improved 
model (1151) is invariant under the first supersymmetry 5^' for arbitrary superpotentials. It 
differs from the naive discretization j8]) of the continuum action (pQ) by the improvement 
term a W x (d(j)) x which turns into an integral over a total derivative in the continuum, 
and hence zero for periodic boundary conditions. The improvement term is needed for an 
invariance of the lattice action under one supersymmetry transformation. For an interacting 
theory the lattice action ( fl~5l ) is not invariant under the other supersymmetry transforma- 
tion with parameter e; an action preserving only this symmetry can be analogously 
constructed, 

^sy = | U<t>? + a E & lr ^- (16) 

x x,y 

It is invariant under the nilpotent supersymmetry transformations 

5 {2) <px = $ x e, 8 {2) i> x = 0, 5 (2 V. = -Le- (17) 

To generate the same fermionic term as in S snsy we choose £ x = — (d(f>) x + W x with antisym- 
metric (d xy ) and symmetric W xy . Then the supersymmetry ( jT7j ) agrees with 5^ in (jHj), and 
the action takes the form 

^SUSy = | ^ {(d(f)) X ~ W x f + a y] j} x {d xy + W xy ) lj)y 

x x,y 

= S naiye -aJ2WM(d<t)) x (18) 
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For periodic fields S snsy and S snsj converge to the same continuum limit. On the lattice they 
are only equal in the noninteracting casein 

2.3.1 The improved models in detail 

We consider three supersymmetric versions of the discretization (fT5l) with improvement 
term. 

(iv) Supersymmetric model with Wilson fermions and Ito prescription 

In order to avoid doublers and at the same time keep half of supersymmetrv we use the 
antisymmetric matrix (d xy ) and shift the superpotential by a Wilson term [30] 13 The corre- 
sponding model 

S$=S snsy , d = d, W x = ~{A<j>) x + W'(<f> x ), W xy = ~A xy + W"(<f) x )5 xy (19) 

is invariant under the supersymmetry 5^ [36j. Of course, the non-interacting model is also 
invariant under S^ 2 ' . 

With these definitions, the improvement term is given by the well-known Ito prescription 

(v) Supersymmetric models with Wilson fermions and Stratonovich prescription 

Instead of the Ito prescription, we can choose the Stratonovich scheme BU] for the evaluation 
of the surface term, ^2 x W'(cr x ) (<fi x — with a x = \{4> x + x -i)Jj The corresponding 

action can be obtained from ([TBI) with a Nicolai variable £ x (<j>) = (d h <f)) x + W'(a x ), 

S[ 5 J = S susy , d=d, W x = -^{A ( j ) ) x + W , {a x ), W xy = -~A xy + dW ^ ax \ (20) 

One should note that this procedure differs from the one proposed in [HHj, where the fermions 
are first integrated out in the continuum theory, and only then a Stratonovich interpretation 
is given for the surface term - in this case, the fermionic path integral of the Euclidean evolu- 
tion operator has to be defined in a non-standard way in order for the bosonic Stratonovich 
Jacobian to cancel the fermion determinant. 



6 In order to preserve the second supersymmetry 5^ also for S susy in l|15p in the absence of interactions, 
its definition will have to be slightly modified only for the Stratonovich prescription to be discussed below. 

7 It is obvious that this is equivalent to working with a shifted lattice derivative as in Q and an unshifted 
superpotential since the action now only depends on the invariant combination £ x . 

8 For monomial superpotentials W[<f) = (f) k , k = 1, 2, . . ., this prescription is equivalent to the prescription 

J2 X \\W'((t>x) + W(<f>x-i)) (<f>x — 4>x-i) ; in the latter case, the superpotential terms are evaluated only at 
a given lattice site. 
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We will see that compared to the fermion determinant involving continuum derivatives, 
the continuum limit of the fermion determinant for the Ito prescription is off by a factor 
depending on the superpotential whereas the Stratonovich prescription reproduces exactly 
the desired continuum result. 

Since fl20l) was constructed from (TT3l) in terms of Nicolai variables, it is manifestly super- 
symmetric under 5^ as given in (TT41) . The discretization of the second supersymmetry is 
in general not preserved on the lattice, not even in the free case. This latter fact suggests a 
modification of £ x (</>) in fTTTl) to 

IM = -(d<f>) x - ^(A0) x + W\a' x ). (21) 

This changes effectively the backward- into a forward-derivative, and the derivative of the 
superpotential is evaluated now at a' x = \{4> x + 4> x +i)- With these definitions, 5^ is a 
symmetry of the action (l20l) in the absence of interactions. It is also this variation with 
which we compute Ward identities in section [3l 

(vi) Super symmetric models with SLAC derivative 

In order to avoid fermion doublers, we can specialize d to be the SLAC derivative, 

S[ 6 J = S Bmy , «9 = «9 slac , W x — W'((f) x ). (22) 

In spite of its nonlocality, the fermion and boson masses extracted from two-point functions 
prove to approach the continuum value quite fast; the quality turns out to be comparable 
to that of the Stratonovich prescription. The interacting supersymmetric model with SLAC 
derivative is only invariant under 5^ in ( |T4l) by construction. 

2.4 Fermion determinants 

In this subsection we demonstrate which sign of the Wilson term we must choose in order 
to guarantee positivity of the fermion determinant. After that, it will become clear that 
as compared to the value for the continuum operator, the fermion determinant for the Ito 
prescription is off by a factor depending on the superpotential whereas the Stratonovich 
prescription and the SLAC derivative reproduce the desired continuum result. 

2.4.1 Sign of the determinants 

For a real fermion matrix d xy + W xy , complex eigenvalues A appear pairwise as A A in the 
determinant. Hence, the determinant can only become negative through real eigenvalues. 
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Since without loss of generality, eigenvectors v x to real eigenvalues A can be taken to be real 
(otherwise, take v x + v*) and normalized, only the symmetric part of the fermion matrix 
contributes to a real A: 

\ = J2 V * ( d *v + W *y) v y = J2 v * ( 9 4 + W *v) v v ( 23 ) 

x,y x,y 

For the antisymmetric SLAC derivative the symmetric part <9 S is absent and no Wilson term 
is required, W xy = W , (<p x )8 xy , such that the real eigenvalues are given by 

A = £V"(&)t& (24) 

X 

We conclude that all real eigenvalues and thus the determinant will be positive for the 
models (Hi) and (vi) in case W" is nonnegative definite. For the models (i), (ii) and (iv) 
with Wilson term the real eigenvalues are given by 

A = -fE"^ + E W "^>1- ( 25 ) 

xy v 

Since —A is positive all real eigenvalues A and therefore the determinant will be positive 
definite in case the W"((p x ) are nonnegative and r > 0; as usual, in this paper we choose 
r = 1 in this case. Vice versa, for a negative W"(<p x ) we would have to change the sign of 
the Wilson term (i.e., choose r = —1) for the fermionic determinant to stay positively In 
the following subsection we shall prove by an explicit calculation that also for the model (v) 
with Stratonovich prescription the fermionic determinant is positive for positive W" . 
If W" is neither positive nor negative definite positivity of the fermion determinant is not 
expected. In fact, for instance in the particular example of W (and W") being an odd power 
of (f) we have to expect a change of sign: From the interpretation of the Witten index 



/ DtpDtpD(j)e~ s = j D(f> det(d + W" ((f))) e~ ShOB . (26) 

J per.b.c. J 

(for the path integral with periodic boundary conditions for all fields) as the winding number 
of the Nicolai map £ = d(p + W regarded as a map from the space of bosonic variables to 
itself [42], we expect that there may be phases with broken supersymmetry for W even; the 
Witten index vanishes. This would be impossible for a determinant with definite sign. 



3 Either way, this part of the action satisfies site- as well as link reflection positivity [H] 
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2.4.2 Calculating the determinants 

The fermionic determinant for the Ito and Stratonovich prescription can be computed ex- 
actly. The (regularized) determinant of the continuum operator d T + W"((/)(t)) on a circle 
of radius (3 can easily be seen to be [13] 



. ( . d T + W"{<j>{r)) \ = ^H\^drW"{<P{r))) _ 



d T + m J sinh(fm) 

this is the value with which we have to compare the lattice results. Note that for a non- 
negative W" the determinant is positive. 

For Ito 's calculus with Wilson fermions, the ratio of the determinant of the fermion matrix 
for the interacting theory 

d xy + W xy = d h xy + W"(<p x )5 xy (28) 

to that of the free theory is given by 

Hpt ( d+(W xy ) \ U^ + aW"^))-! 

aei y d b + ml j ito (1 + am)N _ 1 ■ y^) 

It is positive for positive W"(4>) and this agrees with the results in the previous section. For 
iV = (3 /a — > oo, the product converges ta^_ 

det (3UW^\ »~ ejW^* det (?^m.) (30) 
V<9 b + ml J lto e?™/ 2 \ d T + m J y ' 

since Infill + aW"(<f> t )) = E ln (! + aW"{<j) x )) -> / dxW"(<f> x ). The limit pD differs by 
a field dependent factor from the continuum result. 

For Wilson fermions with Stratonovich prescription the regularized determinant of the 
fermion matrix 

d xy + W xy = d h xy + l -W"{a x ) (5 xy + 5 x . 1>y ) (31) 

is again positive for positive W" . But in contrast to the determinant ( 1291 ) with Ito prescrip- 
tion it converges to the continuum result, 

det U + (m x ,)J strat n(i+»-ri(i-» v d T+m } {62) 

where m xy = \{b~ xy + 6 x -x,y) and M x = W"(a x ). One can show that for N — > oo the 
fermionic determinant with SLAC derivative converges rapidly to the continuum result. 

10 In a completely analogous manner, the right-derivative would lead to the inverse prefactor in front of 
the continuum result. 
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3 Simulations of supersymmetric quantum mechanics 



We have performed high precision Monte-Carlo simulations to investigate the quality of 
the six lattice approximations introduced in subsections 12.2.11 and 12.3.11 The models with 
actions S^J , S^J and S^J are not supersymmetric whereas the actions , S[ 5 J and S^J 
preserve the supersymmetry 



3.1 Effective Masses on the lattice 

In order to determine the masses we have calculated the fermionic and bosonic two point 
functions 

= (4M and Gjl^HW (33) 

in all models in) and fitted their logarithms with a linear function. This way of determining 
the masses mbos(») and mf erm (a) from the slope of the linear fit works well for all models 
with ultralocal derivatives. Details on technical aspects of the extraction of effective masses 
can be found in section 15.31 



3.1.1 Models without interaction 

The actions of the non-interacting models are quadratic in the field variables, 

Sfree = ^ ^ <f> X K xy (f>y + ^ 1p x M xy 1p y . (34) 
xy xy 

Actually, for W'(<f>) = m<j) all the improvement terms vanish and only four of the six actions 
introduced in the last section are different. The corresponding matrices M and K are given 
in the following list0 





°free 


o(4) _ o(2) 
free free 


q(5) 
°free 


<M _ c(3) 
free free 


M 


d b + m 


d b + m 


(1 - ?f)d h + m 


C^slac + m 


K 


-A + m 2 


—A + m 2 — amA 


-(l-(af) 2 )A + m 2 


-(9 8lac ) 2 + m 2 






= M T M 


= M'M 


= M T M 



(35) 



For the actions in the last three columns we have detK = (detM) 2 as required by su- 
persymmetry. This is not true for the non-supersymmetric naive model S^ e with Wilson 
fermions. 

"Note that d h + (d h ) T = -A and -d 2 + ±a 2 A 2 = -A. 
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Without interactions, the masses of all models as determined by numerically inverting M 
and K converge to the same continuum limit m. For the free supersymmetric models (ii)- 
(vi), m hos (a) and mf erm (a) roughly coincide even for finite lattice spacings. This is to be 
expected for supersymmetric theories. In the first model, mb os (o) for finite a is already very 
close to its continuum value, in contrast to mf erm (a) which for finite a is considerably smaller 
than m^ os (a). Since the free supersymmetric action S^ e has the same fermionic mass as 
S^ c and since m bos (a) m fcrm (a) for this model, we conclude that its bosonic mass for finite 
a is notably smaller than its value in the continuum limit. Thus, when supersymmetrizing 
the naive model with Wilson fermions we pay a price: the boson masses get worse while 
approaching the fermion masses. 

The situation is much better for the other supersymmetric models with SLAC derivative 
or Wilson fermions with Stratonovich prescription. The masses are equal and very close to 
the continuum result already for finite a. The masses for the Stratonovich prescription are 
comparable to the boson masses of the naive model without supersymmetry. The masses 
for the free models and their dependence on a are depicted in figure [U 



10 



9.6 



1 9.4 



9.2 



(SI I Si 



3 * 



0.005 



0.01 0.015 

lattice spacing a 



0.02 



l bos 

l W 
"Term 

(5) 
1 ferm 

(6) 
l ferm 



0.025 



Figure 1: The masses determined by numerically inverting the kinetic operators for the free 
theories. There are only 4 different masses, cf. (l35l) . 
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Figure 2: Boson and fermion masses for Wilson fermions without improvement (model (i)) 
and non-supersymmetric SLAC fermions (model (in)). The parameters for the linear fits 
can be found in table ED 



From ( 1331) we conclude that 




The mass rn\J m (a) extracted from G\J m is close to the continuum value m such that 

"£*"(«)=•"£?(«) (37) 

and this simple relation explains why the linear fit through the masses Tn^^^a) marked 
with red dots in figure CD has such a large negative slope. 

3.1.2 Models with interaction 

We have calculated the masses for the interacting models with even superpotential 

W( ( j)) = ^(f) 2 + ^(f) i W'(<j)) 2 = m 2 <f) 2 + 2m<70 4 + g 2 ^. (38) 

Since in the weak coupling regime the results are comparable to those of the free models we 
have simulated the models at strong coupling. In order to compare our results with those 
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of Catterall and Gregory in [36] we have picked their values m = 10 and g = 100 for which 
the dimensionless ratio g/m 2 equals unity. The energy of the lowest excited state has been 
calculated by diagonalizing the Hamiltonian on large lattices with small a and alternatively 
with the shooting method. With both methods we obtain the continuum value 

m phys = 16.865. (39) 

We now summarize the results of our MC simulations. As in the free case, m bos (a) ^ 
^ferm(cO for the non-supersymmetric model with action S^J , see figure [3 In addition, for 
g 7^ their continuum values are different and none of the two values agrees with ([39]) . This 
has been predicted earlier by Giedt et al. [37]. We conclude that the naive lattice model 
with Wilson fermions is not supersymmetric for a — > 0. 

The Monte-Carlo results are much better for the second model with action as given 
in ( ITT! ). Although this model is not supersymmetric, its boson and fermion masses are 

(2) 

almost equal for finite lattice spacings. Linear extrapolations to vanishing a yield ?^ bos (0) = 
16.68 ± 0.05 and rn^ m (0) = 16.73 ± 0.04 which are quite close to the correct value 16.865. 

(2) 

The results for the non-supersymmetric model with action S{ d and the supersymmetric 
model with action are almost identical, similarly as for the free models. The masses 
for various lattice constants between 0.005 and 0.03 for the two models are depicted in 
figure El The corrections to the continuum value are of order 0(a) and are as big as for the 
corresponding free models. The slope and intercepts for the linear fits are listed in table [TJ 
At finite lattice spacing a, the masses mb os ,ferm(fl) for model (v) with Wilson fermions and 
Stratonovich prescription are much closer to their respective continuum limits than for 
the model (iv) with Ito prescription. Furthermore, the extrapolated continuum masses, 
m bos(0) = 16.78 ± 0.04 and Tn^ m (0) = 16.77 ± 0.02, are very close to the correct value 
(|39l) . The data points for the supersymmetric models with Wilson fermions are depicted in 
figure [H Again the slope and intercepts of the linear fits can be found in table [TJ Of all 
lattice models with ultralocal derivatives considered in this paper this model yields the best 
predictions. 

The model with Stratonovich prescription for the improvement term is outperformed only 
by the models (Hi) and (vi) with nonlocal SLAC derivative for fermions and bosons. This 
observation is not surprising, since the remarkably high numerical precision of supersymmet- 
ric lattice models with SLAC derivative has been demonstrated earlier in the Hamiltonian 
approach in [9]. Furthermore, this is in line with our results for the free models, see figure [TJ 
The masses for the interacting unimproved model (in) are plotted in figure [2] and those 
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Figure 3: The masses mb os (a) and mf erm (a) for models (it) (shifted superpotential) and (iv) 
(Ito improvement). 
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Figure 4: Masses for supersymmetric models (vi) (Ito prescription) and (v) (Stratonovich 
prescription). Only lattices with at least 61 sites are included in the fit. 
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Figure 5: Masses for the supersymmetric model (vi) with SLAC derivative as compared to 
the masses of the model (vi) with Wilson fermions and Ito prescription. 



for the improved model (vi) in figure El Even for moderate lattice spacings the masses 
m bos( a ) ~ m ievm( a ) are very close to their continuum limits m^ s (0) = 16.84 ± 0.03 and 
m ferm(0) = 16.81 ± 0.01 which in turn are off the true value 16.865 by only some tenth 
of a percent. The extrapolated masses for the unimproved lattice model in table [T] have a 
comparable precision. But of course there is no free lunch, since for the SLAC derivative 
one must smooth the two-point functions Gbos.ferm^) with a suitable filter for a sensible 
mass extraction. Details on the filtering can be found in subsection 15.31 
In the following table we list the slopes k and intercepts m(0) of the linear fits 

m hos (a) = &bos ■ a + m bos (0) and m ferm (a) = A; fcrm • a + m ferm (0) (40) 

to the measured masses for the six lattice models considered. 

The linear fits for the models with improvement (iii)-(vi) are compared in figure [6l Lattice 
supersymmetry guarantees that the boson and fermion masses are equal for these models. 
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model 


^bos ^bos(O) 


^ferm ^ferm(O) 


<j(i) 

o(2) 
D ld 

o(3) 


139.52 ±8.45 12.23 ± 0.08 
-136.85 ± 5.22 16.68 ± 0.05 
-25.22 ±6.24 16.92 ± 0.07 


-186.25 ±4.98 18.40 ± 0.05 
-146.10 ±3.84 16.73 ±0.04 
-33.64 ±2.52 16.97 ±0.03 


c(4) 
°ld 
c(5) 
D ld 
0(6) 


-135.11 ± 7.36 16.68 ±0.07 
-40.40 ±4.46 16.78 ±0.04 
-17.97 ±2.41 16.84 ±0.03 


-138.50 ±2.85 16.64 ±0.03 
-37.55 ± 1.98 16.77 ±0.02 
-18.53 ±0.91 16.81 ±0.01 



Table 1: Slope and intercepts of linear interpolations for the masses. 




Figure 6: Linear fits to the masses of all three supersymmetric lattice models. 



3.2 Ward identities 

The invariance of the path integral measure under supersymmetry transformations leads to 
a set of Ward identities connecting bosonic and fermionic correlation functions. Namely, 
the generating functional for Green's functions should be invariant under supersymmetry 
variations of the fields, 

= 5Z[J,9,9} = I P(0,VOe _5+E * (J ^ +e ^ + ^^ (41) 
J y 
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Figure 7: Ward identities for the free theory, model (ii) (with shifted superpotential). 

Ward identities are obtained from derivatives of this equation with respect to the sources. 
The second derivative d 2 /d.J x dO y leads to the Schwinger-Dyson equation 

(^ y 6^S) = (<M (1) Vv) + <<M (1 V*>, (42) 

whereas d 2 /dJ x d9 y yields 

(M y 5^S) = (<M( 2 V y > + (V^ (2) <M- (43) 

For the improved models (iv)-(vi) in subsection 12.3.11 the actions are manifestly invariant 
under S^ 1 ' , and the left-hand side of (T4*2l) vanishes. Thus for these models the following Ward 
identities hold on the lattice: 

(My) ~ (My) = 0. (44) 

This can be confirmed in numerical checks and merely serves as a test bed for the precision of 
the algorithms. The discretization of the second continuum supersymmetry transformation, 
however, only leaves these lattice actions invariant in the free case. With interactions, the 
term 5^S leads to a nonvanishing left-hand side in the Schwinger-Dyson identities (l43l ) 
which therefore measures the amount by which the second continuum supersymmetry is 
broken by the discretization. Since this supersymmetry can be made manifest by choosing 
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Figure 8: Ward identities for the theory at strong coupling, model (%%) (with shifted super- 
potential). 



different Nicolai variables £, the supersymmetry breaking terms are the difference of both 
actions (TT5l) and (TT8l) . i.e., the discretization of the surface terms (for (iv) and (vi)). 
For the actions in (i)-(iii) all supersymmetries are broken by the discretization; the cor- 
responding Schwinger-Dyson equations again can serve as a measure for the quality of the 
lattice approximation to supersymmetry. Barring interactions, we expect them to hold for 
both supersymmetry transformations in (it) and (in). 

3.2.1 Ward identities of the unimproved models 

The Ward identities have been simulated for models with the same superpotential as in (l38l) . 
Since at weak coupling similar results can be expected as in the free theory with g = 0, we 
have also studied the theory at strong coupling with g = 800 and m = 10 corresponding to 
a dimensionless ratio g/m 2 = 8. 

The naive discretization (i) is not supersymmetric even without interactions, therefore we 
concentrate on the other models where the supersymmetry of the free theory sets a scale for 
the quality of the simulation of the broken supersymmetries in the interacting case. For the 
unimproved models, we can use translation invariance and measure the right-hand sides of 
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Figure 9: Ward identities for the free theory, model (iv) (Wilson fermions with the Ito 
prescription). 

the Ward identities ([42]) and g3J(, i.e. 

n£l y = ty x $ y ) - (Md<P)y) - (<f> X Wy) (45) 

and 

KSly = (<f> X (d<f>) y ) - (0 X Wy) ~ (^y) (46) 

as functions of x — y. It should be noted that without interactions, the first Ward identity 
reduces to a matrix identity D~ l — (D T D)^ 1 D T = for the free Dirac operator D xy = d xy + 
mS xy if one uses that (ip x ipy) = D~ y l and (<f> x <t>y) — {D T D)~y. The corresponding data for 
the free theory in the case of the model with shifted superpotential (model (ii)) are shown in 
figured In order to keep statistical errors small, in all situations 4 runs with 10 6 independent 
configurations were evaluated. Within our numerical precision, supersymmetry is broken for 
this model if the simulation data in the interacting case exceeds the bounds set by the free 
theory. The results for (1451) and ( 1461) at g = 800 are displayed in figure El Remarkably, the 
statistical error is much smaller than in the free theory; supersymmetry breaking is equally 
strong for both supersymmetries S^ 1 ' and 5^ from ([91) at strong coupling. 
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3.2.2 Ward identities of the improved models 

For Wilson fermions with exact supersymmetry and the Ito prescription (model (iv)), the 
statistical error as measured by the free Ward identities is roughly of the same size as for the 
unimproved model (ii), cf. figure El Again, this was obtained by 4 runs with 10 6 independent 
configurations. As expected, the first supersymmetry 5^ (cf. ([HI) and figure [Toll is preserved 
even at strong coupling; however, the supersymmetry breaking effects for 6^ are about three 
times as large as for the corresponding svmmetrv in the unimproved situation. 

o Ward-Identity 1: m = 10, g = 800 
Ward-Identity 2: m = 10, g = 800 

* . _^ 

3jE 

IE 

UUUU 5 10 15 20 

lattice point x 

Figure 10: Ward identities for the theory at strong coupling, model (iv) (Wilson fermions 
with the Ito prescription). 



0.004 



0.002 



-0.002 



-0.004 



For the free supersymmetric model with Stratonovich prescription (model (iv)), the second 
supersymmetry (IT7I) is only a symmetry with the definition (|2"T1) of £. The Ward identities 
of both supersymmetry transformations, 



- (MWv) - (fawx****^)) = o, 

(Mty)v)-(&W'(*e±£*±))-@ x 1> v ) = 0, (47) 
reduce in the free theory with W'((j)) = mip to matrix identities for the free Dirac operator 

D X y 9 X y ~\~ ~2\3%y ^x,y—\)i 

D~ l - (D T D)- 1 D T = 0, 
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Figure 11: Ward identities for the free theory, model (v) (Wilson fermions with the 
Stratonovitch prescription). 
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Figure 12: Ward identities for the theory at strong coupling, model (v) (Wilson fermions 
with the Stratonovitch prescription). 
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-{D T D)- 1 D + (D~ 



0. 



(48) 



Here, the second identity holds since D is circulant and therefore normal. The corresponding 
plot of the left-hand sides is shown in figure HH This determines the mean error above 
which we take supersymmetry to be broken if we switch on interactions. In figure [El 
the first supersymmetry is preserved within a high numerical accuracy whereas the second 
supersymmetry is clearly broken by effects about three times the size of the supersymmetry 
violation in the model with Ito prescription. 

For the supersymmetric model with SLAC derivative, the Ward identities are satisfied within 
statistical error bounds, but determining the mean error in an analogous manner fails since 
the contributions of the bosonic two-point functions ((j) x (d<p) y ) in the free theory lead to 
large errors which obscure the interpretation of the corresponding Ward identities. At 
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Figure 13: Ward identities for the theory at strong coupling, model (vi) (supersymmetric 
model with SLAC fermions). 



strong coupling, however, the interaction terms dominate the action, and the fluctuations 
of the bosonic propagators become increasingly less important so that we obtain rather 
precise results at g = 800. Within the error bounds of model (v), the first supersymmetry 
is obviously preserved in the interacting theory (cf. fig. [[3D, whereas the supersymmetry 
breaking effects of the second supersymmetry are about as large as for Wilson fermions with 
the Stratonovich prescription. 
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4 The J\f = 2 Wess-Zumino model in 2 dimensions 



In this section, we study different discretizations of the two-dimensional Wess-Zumino model 
with (2, 2) supersymmetry. The minimal variant contains a Dirac spinor field and two real 
scalar fields ip a which are combined into a complex scalar field = (fx + if 2 . Also, we use 
the complex coordinate z = x l +ix 2 and its complex conjugate z in Euclidean spacetime and 
denote the corresponding derivatives by <9 = \{d\ —182) and <9, respectively. The Euclidean 
action contains the first and second derivatives W and W" of a holomorphic superpotential 
W(4>) with respect to the complex field 0, 

S cont = J d 2 x (2<90<90 + ^ I W'\ 2 + ipMi?) , M = + W"P+ + W"P„ , (49) 

where P± = |(1 ± 73) are the chiral projectors. In the Weyl basis with 7 1 = a\, 7 2 = — a 2 
and 73 = «7 X 7 2 = (T3, the complex spinors can be decomposed according to 

) and ^ = (-01,-02). (50) 



In this basis, the supersymmetry transformations leaving S^nt invariant are 

50 = fret + 6ij l = -\W'e x - 90e 2 , Si/; 1 = -\W'e x + d^e 2 , 

50 = ^e 2 + e 2 tfj 2 , dtfj 2 = -<90£i - \W'e 2 , Stfj 2 = - \We 2 . l ° 

Similarly as in quantum mechanics a naive discretization of this model breaks all four super- 
symmetries. In order to keep one supersymmetry one can add an improvement term. In what 
follows we shall only consider improved models; they differ by our choice of the lattice deriva- 
tives. Instead of trying to generalize the Stratonovich prescription to the two-dimensional 
situation, we find in subsection 14.11 that a non-standard choice of the Wilson term leads to 
an improved behavior in the limit of vanishing lattice spacing. This is corroborated by the 
results of our simulations for the case of a cubic superpotential W = |m0 2 + |g0 3 which 
we present in subsection 14.31 

4.1 Lattice models with improvement 

We start with a two-dimensional periodic iVi x N 2 lattice A with complex bosonic variables 
<p x and two complex spinors ip x ,ip x on the lattice sites x = (x 1 , x 2 ) G A. Again two different 
antisymmetric lattice derivatives in direction /x are used. These are the ultralocal derivative 

9/j,,xy = ^ (<5x+e M ,j/ — ^x-e M ,y) (52) 
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with doublers and the nonlocal SLAC derivative without doublers, which for odd Ni,N 2 
reads 



aslac _ oslac r oslac _ oslac r , oslac _ n /ro\ 

u l,xy — u x 1 ^y 1 °X2,y2) u 2,xy ~ u x 2 ^y2° x i^ dI1U U fJ.,xx ~ u V' JO / 

(analogously to the one-dimensional SLAC derivative defined in (jlj)). Later we shall remove 
the fermionic doublers of by introducing two types of Wilson terms, both containing 
the lattice Laplacian A defined by 

(A0)(x) = 4( E [0(^ + e M ) + 0(x-e M )]-40(x)). (54) 

As for the continuum model we use the holomorphic lattice derivative d xy := \{di jXy — id^ xy ). 
The Nicolai variables of one-dimensional systems are easily generalized to two dimensions, 

& = 2(90), + W*, & = 2(90), + W^; (55) 

again, W x denotes terms (to be specified below) derived from the superpotential. The 
bosonic part of the action is Gaussian in these variables, Sb OS = ^-z2 x ^x^x, an( ^ nas the 
explicit form 

S hos = a 2 J2 (2(d$) x (d<f>) x + W x (d<P) x + W X (B^) X + \\W X \ 2 ) . (56) 

X 

For antisymmetric derivatives d^, the kinetic term has the standard form ^ Q „_ 1 (9 /Lt </? )^ in 
terms of the real fields (p a ; in particular this holds true for the ultralocal derivative <9 M and 
the SLAC derivative <9^ lac introduced above. The second and third term in S^ os are absent in 
the continuum action. In a naive discretization of the continuum model these improvement 
terms do not show up. In the continuum limit they become surface terms and could be 
dropped. On the lattice they are needed to keep one of the four supersymmetries intact. 
Supersymmetry requires an additional fermionic term Sf erm = a 2 ^2 ip x M xy i/j y for a two- 
component Dirac spinor field in such a way that the determinant of the Jacobian matrix 

dC/tyy <K,:<)oA (W xy 2d xy \ = dw x 

dUd<f> y d&/d$ v J \2d xy ir,J ' xy <k hl ' 1 

for the change of bosonic variables (0,0) i— > (£,£) cancels the fermion determinant detM. 
Actually the fermionic operator M in (T49l) with 7-matrices in the Weyl basis and continuum 
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derivatives replaced by lattice derivatives is identical to the Jacobian matrix. Hence we 
choose as the fermionic part of the action 

Sferm = 0? ^ ^M^y, M xy = Y^xy + W xy P + + W xy P_ 

= M + W"((f) x )8 xy P + + W"{(j) x )8 xv P_ . (58) 

By construction the action S snsy = Sb os + 5ferm with improvement terms is invariant under 
the supersymmetry transformations generated by 

8w<t> x = m>i, s^m = -i^, s^i = 0, 



this corresponds to a discretization of the continuum symmetry (1511 1 with E\ = e 2 = and 
E\ = 82 = s. The other three continuum supersymmetries are broken; for appropriately 
chosen Nicolai variables £ an actionpreserving any one of the three other supersymmetries 
can be constructed analogously [30jo In this paper, we are going to use the Nicolai variable 
(l55l) and consider several possibilities to remove fermion doublers. 



4.2 The lattice models in detail 

We introduce three different lattice approximations to the continuum Wess-Zumino model 
( 1491) . They are all equipped with an improvement term and thus admit one supersymmetry. 
The first two models contain Wilson fermions and the third the SLAC derivative. It will 
turn out that the discretization errors of the eigenvalues of the bosonic and fermionic kinetic 
operators in the free case indicate how good the approximation to the continuum theory is 
when we turn on interactions. 

(i) Super symmetric model with standard Wilson term 

Here we choose ultralocal derivative <9 M and add a standard Wilson term to the superpotential 
to get rid of the doublers of 7^(9^ so that 

= S b os + 5 fcrm with d, = d„ W x = -^-(A<p) x + W'(<p x ). (60) 

For later convenience we do not fix the Wilson parameter r in this section. In this lattice 
model the Dirac operator M in (l58l takes the form 

= Yd, ~ yA (61) 

12 The corresponding Nicolai variables can be read off from the right-hand sides of 5ip a in ([5T]) for e\ = 
±S2 = e and e = or from the right-hand sides of 5"tp a for E\ = ±e 2 = e and e = 0. 
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and we easily recognize the standard Dirac operator for Wilson fermions. The bosonic part 
of the action may be expanded as 

S£L = J E (M K ^ + \W'M\ 2 ) +a 2 E (W'M (80 - ^A0) x + c.c.) (62) 

with the kinetic operator 

K = -A + (±arA) 2 , where A = -<9 M d M . (63) 



The last term in (|62l) is the improvement term — a discretization of a surface term in 
the continuum theory. Note that even for the free massive model with W'((j)) = mcj) the 
improvement term — |(amr) a 2 ^ (px(A(j)) x is non-zero such that (l62l becomes 

t 
2 



S bi = \ Yl ^ K xy<l>y with kW =K + m 2 - armA. (64) 



xy 

^2 „ J ^2 



The eigenvalues of the commuting operators —A and —A are p and p , where 

1 2 fO'PuX 2k ku 

p^ = -smap^ Pn = - sm ^— J with p M = — fc M G {1, 2, . . . , JV M }, (65) 

A* 

such that the eigenvalues of the matrix IfW in (1641 are 

H P = p 2 + {m + \arp 2 ) 2 . (66) 

On the other hand, the eigenvalues of the free Dirac operator Mq 1 ^ + m are given by 

Ap = m + |arp 2 ± i\p\. (67) 

Thus, /i p = Xt X~ and the fermionic and bosonic determinants coincide for the free theory. 
The discretization errors for these eigenvalues are of order a, 

jj p = p 2 + m 2 + (arm)p 2 + 0(a 2 ), Aj = ±i\p\ + m + \ar p 2 + 0(a 2 ). (68) 

It is remarkable that the bosonic part of the action in the continuum is an even function of 
the mass, whereas its discretization (1621 ) is not. This spoils the sign freedom in the fermion 
mass term on the lattice and motivates the following second possibility. 

(ii) Super symmetric model with non-standard Wilson term 

Again we choose the ultralocal derivatives <9 M but now add a non-standard Wilson term to 
the superpotential to get rid of the doublers, so that the action now is 

S^ = S hos + S ierm with d, = d„ W x = % -^(A<P) X + W\<p x ). (69) 
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This choice is equivalent to an imaginary value of the Wilson parameter inside the holomor- 
phic superpotential. Now the Dirac operator M in (l58l) has the form 

M (2) = + ^73A (70) 

and contains a non-standard Wilson term, reminiscent of a momentum dependent twisted 
mass. It should be noted that the only difference between the bosonic actions (1621 ) and 

S£ = y E + \W'M\ 2 ) +a 2 J2 (Wfa) (d<j> - ^A0) x + c.c.), (71) 

X X 

is the improvement term. The modified Wilson term in f l69l) yields an action which is even 
in the mass m. Actually, for the free massive model with W'(<f>) = mcf) the improvement 
term vanishes and 

4I = yE^ with K^ = K + m\ (72) 

xy 

The eigenvalues of and of the free Dirac operator + m in this case are given by 



fx p = m 



+ p 2 +(^arp 2 ) 2 and A^ = m ± i\J p 2 + Qarp 2 ) 2 . (73) 



Again the determinants of fermionic and bosonic operators are equal. The added advantage 
in this situation is, however, that the discretization errors of the continuum eigenvalues are 
only of order 0(a 2 ). Namely, for small lattice spacing, 

/i p = m 2 +p 2 + ko 2 + 0(a 4 ), Aj = m ± i(p 2 + Ka 2 ) 1/2 + 0(a A ), (74) 

where the 0(a 2 )-coefficient 

K = + j(£4 2 = (j - 1) Erf + Tif^f^) 2 < 76 > 

vanishes for 3r 2 = 4 and pi = or p 2 = 0. As explained in section [5731 we take spatial aver- 
ages over two-point functions for the mass extraction; this projects the spatial momentum 
to zero. Thus, one might expect discretization errors of order 0(a 4 ) for 

r 2 = t. (76) 

In fact we will see in section 14.31 that this choice leads to the best continuum approximation 
- and this in spite of the fact that it violates reflection positivity (as does the improvement 
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term in all supersymmetric models). 

(Hi) Supersymmetric model with SLAC derivative 

The Dirac operator ^cH ac with nonlocal and antisymmetric SLAC derivatives defined 
in f l53l) has no doublers and no Wilson terms are required; this leads to the action 

= S hos + S icrm with d, = d^ c , W X = W'(4> X ) (77) 

with Dirac operator 

M (3) = 7 ^^ lac , (78) 

cf. (l58lh For the free massive model with SLAC derivative the improvement term vanishes, 
in particular the bosonic part of is 

fi = yE^* with ^ (3) = -A slac + m 2 ; (79) 

xy 

all supersymmetries are realized. We shall see in section 14.31 that is a very good ap- 
proximation to the continuum model and in section [6] that the lattice model based on the 
SLAC derivative is one-loop renormalizable in spite of its nonlocality. 

4.3 Simulations of the Wess-Zumino model 

In subsection 14.11 we have formulated the model in a complex basis, which is natural and 
convenient for models with two supersymmetries (in particular, the simplest form of the 
Nicolai map (1551) is in terms of the complex scalar fields <fi — (fx + iy?2 and £ = £i + i^). 
On the other hand, for numerical simulations it is convenient to have a formulation of the 
model in terms of the real components (p a and £ a which are combined to real doublets, 

V=( lfl ) and H=(i 1 ). (80) 



\<P2/ V 6 / 

As to the fermions, it is most appropriate to use a Majorana representation with real 7- 
matrices 7 1 = cr 3 , 7 2 = o~\ such that 73 = ry x 7 2 = — o~ 2 - All simulations were done for the 
model with cubic superpotential W = ^mcj) 2 + |g0 3 with derivative 

W'((p) — m<fi + u + i(m<f2 + v), (81) 

where we have introduced the abbreviations u = g{^\ — fl) and v = 2g{p 1 <p 2 - The actions 
contain a quartic potential 

V(i P ) = (g 2 <p 2 + 2mg<p 1 )p 2 . (82) 
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In terms of real fields, the Nicolai map (1551) takes the form 

£(«) = (M (n) + m)cp + ^ for n = 1,2,3 (83) 

with model-dependent free massless Dirac operators Mq as given in (IBTl ). ( 1701 ) and f78| . 
The bosonic actions for the models with standard Wilson, modified Wilson and SLAC 
fermions can now be written as 

St = y K SV.) + y £ ^vO + AW ( 84 ) 

with model-dependent kinetic operators lO n ' as introduced in (1641) . (1721) and (|79j) . In terms 
of the divergence and curl of a vector field in two dimensions, div ip = d\(f \ + d 2 <f2 and 
curl (p = di(f2 — d2<fi, the model-dependent improvement terms 



A(1) = fl2 £ u * (( div ^ ~ y( A< ^i)-) ~ fl2 £ ^ (j curl v)* + y ( A ^)^ 

X X 

A< 2 > = a 2 E^((div(^) x -y(A^ 2 ),) +a 2 E^(-(curl^) :c + y(A^ 1 ) 1 .), (85) 

X X 

A (3) = a 2 ^ ^x(div y)x - a 2 £ ^x(curl cp) x 

X X 

again are discretizations of continuum surface terms. 
4.3.1 Models without interaction 

As expected, the masses of all models without interactions converge to the same continuum 
limit m. Since all free models are supersymmetric (w.r. t. two supersymmetries), boson 
and fermion masses extracted from the two-point functions coincide even at finite lattice 
spacing. The masses m(a) for the models (it) and (Hi) with non-standard Wilson term 
(with r 2 = 4/3) and the SLAC derivative, respectively, at finite a are already very close to 
their continuum limits, the effective mass as a function of a for the model with standard 
Wilson term (with r = 1) has a much larger sloped This is in line with the approximation 
of the eigenvalues (l68l) and (1741) to those of the continuum kinetic operators. 



13 This behavior is reminiscent of the masses m(a) for the corresponding quantum mechanical model in 
section l|3.1.ip . 
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4.3.2 Models with interaction 



We have calculated the masses for the interacting models with cubic superpotential 



for masses m — 10 and couplings g ranging between and 1. The effective masses at 
different values of the lattice spacing are determined as described in section I5T31 Again, due 
to supersymmetry boson and fermion masses coincide also for finite lattice spacing. In all 
cases, they converge to a continuum value which however cannot reliably be distinguished 
from the value mf ree (a = 0) = 10 within error bounds. This is to be expected since in 
continuum perturbation theory, the one-loop corrected mass is 



in a renormalization scheme without wave-function renormalization (this corresponds to a 
correction less than about 0.3% with our values of m and g). Thus, significant effects should 
only be seen at larger values of g/m. Unfortunately, our simulations require reweightings 
which lead to rather large error bounds (which grow as the coupling increases). The conver- 
gence behavior to the expected continuum value is model-dependent: The model with the 
non-standard Wilson term shows the expected improved behavior leading to good estimates 
for the continuum mass already at finite lattice sizes. As in quantum mechanics, this applies 
also to the SLAC derivative. 



5 Algorithmic aspects 

In this section we briefly outline the methods and algorithms we have used in our simulations. 
In particular we have modified the treatment of the fermion determinant in the well known 
Hybrid Monte Carlo algorithm (HMC) [44] . 

Although due to the low dimensionality of the models under consideration the numerical 
studies to be carried out are much less demanding than, e. g., four-dimensional lattice QCD 
we have to face similar problems with respect to the treatment of the fermion fields. However, 
the computational tasks at hand allow for strategies which are more accurate and easier to 
implement than what is widely used there. Since our models do not contain any gauge 



W{4>) = -m0 2 + -g4> 3 



(86) 




(87) 
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Lattice spacing 


Wilson 


Twisted Wilson 


m B 


mp 


m B 


mp 


0.1250 


6.34(4) 


6.4872(1) 


9.51(8) 


9.95(7) 


0.0833 


7.32(3) 


7.2730(2) 


10.0(1) 


10.0182(1) 


0.0625 


7.81(8) 


7.768(1) 


10.4(1) 


9.99(1) 


0.0500 


8.16(4) 


8.07(1) 


9.82(1) 


9.93(3) 



Lattice spacing 


SLAC 


m B 


mp 


0.0769 


9.9(2) 


10.0(2) 


0.0667 


10.0(1) 


10.0(1) 


0.0526 


10.0(1) 


10.00(5) 


0.0400 


9.95(6) 


9.99(2) 


0.0323 


10.03(3) 


9.98(1) 


0.0213 


9.83(3) 


9.98(1) 



Table 2: Comparison of extracted masses at g = 0.5 
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standard Wilson boson 
standard Wilson fermions 
SLAC boson 
SLAC fermion 
non-standard Wilson boson 
non-standard Wilson fermion 
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Figure 14: Effective mass for the two-dimensional Wess-Zumino model at g = 0.5 
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degrees of freedom, the resulting Dirac operators have a rather simple form, a fact which 
can be put to practical use as will be shown below. 

Nonetheless, the mere presence of the fermion determinant in the partition function intro- 
duces a nonlocality which has to be taken into account in the Monte Carlo algorithm. We 
thus decided to base our numerical simulations on the HMC which allows for a simultane- 
ous update of all bosonic field variables. For the quantum mechanical models discussed in 
section [2] substantial improvements can be achieved if it is possible to compute the fermion 
determinant in closed form. For the two-dimensional models, however, comparable results 
are not available, and these theories are plagued by a strongly fluctuating fermion determi- 
nant. Even worse, the fermion determinant may take on positive and negative values which 
drives the simulations directly into the so-called sign problem. To proceed we will therefore 
treat the quantum mechanical models again separately from the two-dimensional models 
and discuss them one at a time. 



5.1 Quantum Mechanics 

Our setup for the bosonic degrees of freedom for the HMC does not differ from the standard 
procedure and is formulated on an enlarged phase space involving the real bosonic field X 
and an additional conjugate momentum field ir x . As usual these fields are propagated along 
the molecular dynamics trajectory by Hamilton's equations 

■ dH dH 

<Px = TTx = -7T— , (88) 

dir x d<p x 

where 

H=\^l + S ^\- ( 89 ) 

X 

Since on the lattice the fermions are already integrated out the expression to be used in 
eq. (T891) is given by 

S[<j)] = S B [<f>] + lndet M [</>]. (90) 

In the standard approach one would now introduce a pseudofermion field \ to obtain a 
stochastic estimate for detM[0] which however will necessarily introduce additional noise 
to later measurements. Hence it would be clearly favorable to take the fermion determinant 
exactly into account. While a direct computation of the fermion determinant at each step 
of a trajectory is also feasible in these one-dimensional theories one can do even better due 
to the simple structure of the Dirac operator. Despite their differences in the bosonic part 
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of the action the models involving an antisymmetric derivative matrix and a Wilson term, 
namely the models (i), (ii), and (iv), share the same fermion Matrix M w [0]. For these cases 
a closed expression for the determinant with Wilson parameter r = 1 is given by (l29l) . 

det M w [cj)) = ]J(l + m + 3g(f>l) - 1. (91) 

X 

It should be noted that with this 'diagonal' nonlocal form of the determinant even local 
algorithms could have been used. A similar expression for the Stratonovich prescription as 
given in (l32l ) also allows for a quick computation of the fermion determinant in model (v). 
For the SLAC derivative, it is easier to again exploit the local structure of the interaction 
terms in the fermion matrix by varying 

5(tr\nM[(j)}) =tr(5M[(f)]M- 1 [(f)}) (92) 

in order to determine the contribution of lndetM[0] = trlnM[0] to the equations of mo- 
tion (l88l) . Since we only consider Yukawa couplings, the variation <5M[0] enters (l92l) as 

= 3 #AA*- (93) 

This holds true if the derivative of the superpotential appears in M only on the diagonal; 
this means that only a single matrix element of M _1 [0] has to be computed for each site 
<Az 3 Obviously, we cannot avoid an inversion of M[(j)\ in order to update all sites. 
Finally let us briefly mention some details of our simulation runs. Since the physical value 
of m as well as the physical volume of the box L were kept fixed at L = 10 • rrT x = 
Na, the lattice spacing a was varied with the number of points N in the lattice. The 
lattice sizes we have considered range from N = 15 to N = 243 corresponding to a = 
0.06 . . . 0.004. We have checked that the measurements were insensitive to finite size effects. 
For all lattice sizes we generated between 250 000 and 400 000 independent configurations 
for the mass extraction; in order to improve the signal to noise ratio for the Ward identities, 
10 6 independent configurations were used (with 4 independent runs in order to minimize 
statistical errors). 



14 If this method was to be applied to the Stratonovich prescription (where we can alternatively use |32|) ). 
two matrix elements would have to be computed since the interaction is off-diagonal by one unit in the 
fermion matrix. 
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5.2 The two-dimensional Wess-Zumino model 

For the two-dimensional models, the fermion determinant is not positive definite, in contrast 
to the situation in one dimension. Indeed a numerical experiment shows that the determi- 
nant has large fluctuations at strong coupling and changes its sign if g/m is larger than a 
certain threshold value of order 0(1). This value depends on the lattice spacing (and on 
m). A closer look at the spectrum of the fermion matrix reveals a second related problem, 
namely the existence of very small eigenvalues. They increase the condition number of the 
fermion matrix by several orders of magnitude and prevent a straightforward application 
of the pseudofermion method. In addition, the bosonic potential |jy'(0)| 2 possesses two 
separate minima; this might lead to complications with respect to ergodicity. In order to 
check at which values of the coupling the standard HMC breaks down, we have performed 
simulations in the weakly coupled regime where the aforementioned contributions from the 
fermionic fluctuations can be taken into account via reweighting quenched ensembles. An 
added advantage is that with this approach, field generations can be generated very quickly; 
on the other hand, larger ensembles are used since reweighting requires higher statistics. 

2.5 



2 



1.5 

g 
a, 

1 

0.5 





-25 -20 -15 -10 -5 5 

R 

Figure 15: Probability distribution of R = ln(detM/ detM ). The stronger pronounced 
the peak the better the statistical errors are under control. The failure of the reweighting 
technique is visible for g > 1. The plotted data was generated from 20 000 configurations 
for the model (in) on a 31 x 31 square lattice. 
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Figure 16: Comparison of the bosonic two-point function Gbos(^) for the quenched and full 
theory at g = 0.5 for the model with standard Wilson term (i) on a 32 x 16 lattice. Similar 
results hold for the other models as well. 

Trial runs indicate that for the volumes considered the reweighting method should work for 
g < 1. Figure [T5l shows that for these moderate couplings the fluctuations of the fermion 
determinant are within two orders of magnitude. For g = 2 the fluctuations span more than 
ten orders of magnitude and the number of relevant configurations becomes ridiculously 
small. 

As pointed out in section [4.3.21 for g = 0.5 the perturbation of the mass is about 0.3%, an 
effect which is clearly not visible in our simulations. However, we can compare quenched 
with reweighted expectation values and check whether they are sensitive to the inclusion 
of dynamical fermions at all. We have found that even in this weakly coupled regime 
the effective masses of bosonic and fermionic superpartners coincide only if the fermion 
determinant is properly taken into account. Otherwise, the two-point functions (and hence 
the masses) deviate considerably from their values in the full theory, as can bee seen in 
figure [161 Clearly, in order to simulate at larger couplings an improved algorithm is needed. 
For example, a comparison with the results in [30J would require a ratio of g/m ~ 0.3, which 
obviously is not feasible with the reweighting method used in this paper. 
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5.3 Measurements and determination of masses 



In this subsection we expand on the extraction of effective masses in our quantum mechanical 
models. They are determined by the asymptotic behavior of the (bosonic or fermionic) two- 
point functions G(x) ~ J2i Cie~^ Ei ~ E °^ x : For x far away from the midpoint N/2 of the 
lattice, the asymptotic expression for the bosonic two-point function is dominated by a 
single mass state with E\ — Eq — m c g, and as a symmetric function it is proportional to 
cosh(m e g(x — N/2)); the fermionic two-point function as a superposition of symmetric and 
antisymmetric parts (w.r. t. the midpoint) [30] is for small x proportional to e~ mcfiX . Thus, 
the masses can be extracted by exponential fits from the corresponding simulation data 



The x-region for the fit should be chosen in such a way that 

(i) the contribution of higher energy states with Ei > E% in the asymptotic expression 
G(x) = J2i Cie~^ Ei ~ E °^ x is negligible, 

(ii) the effect of the improvement terms, which violate reflection positivity (i. e., the q are 
not necessarily positive) and therefore damp the two-point function for small values 
of x below the continuum value, better not influence the result, 

(iii) the errors (which are larger for smaller values of the two-point function) should be 



(iv) and in the case of bosons, the influence from the second exponential tail in cosh(m e ff(x— 



N/2)) = i( e m -«( x - N /^ + e -™eff(*-A72)) s hould not interfere with the exponential fit to 



Thus, we are supposed to choose a region for x in between the left boundary (x ^> by 
(i) and (ii)) and the midpoint (x <C N/2 by (iii) and (iv)). This works reasonably well 
for Wilson fermions (the results are presented in section l3?fl) . but for SLAC fermions, we 
observe an oscillating behavior of the two-point functions near the boundaries which requires 
a more careful investigation. 

In order to understand the nature of this phenomenon, we compute the propagator of a free 
massive fermion in the continuum, 





(94) 



minimized, 



the first. 




(95) 



A 
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where the sum runs over the nonzero part of the spectrum of the differential operator <9 + m, 
and ip\ are the eigenfunctions of d + m with eigenvalue A. If we impose periodic boundary 
conditions with ip\(x) = ip\(x + N), we have ip\(x) = -^h*e^r kx with A = + m for all 
integer k. Inserting this into (l95l) . we obtain 

for the two-point function on the circle. If we discretize the circle (and so introduce a 
momentum cutoff), the sum in ( 1961 ) is truncated to a finite number of terms and reduces 
just to the propagator of the SLAC derivative. This cutoff in momentum space leads to the 
Gibbs phenomenon which is in fact what we observe, e. g, in figure [TZl 



0.9 r 




lattice point x 



Figure 17: The free fermionic two-point function (g p hy S = 0, m p h ys = 10) before and after 
the application of the filter. 

The error made by an approximation to a continuous and periodic function by n of its 
Fourier modes decreases exponentially with n. For a discontinuous function, the error is 
proportional to some power n~ 5 away from the discontinuities. This can be improved by a 
filter which increases the convergence rate 5 or even recovers the exponential approximation 
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Figure 18: This plot shows one example of a fermionic two-point function of the interacting 
theory (m p h ys = 10, g p h y s = 100) in a logarithmic representation before and after the 
application of the filter. A linear fit of the filtered function yields the effective mass m. In 
the last part of the graph a large numerical error can be observed. 



of a continuous function^ For our computations, we have used the "optimal filter" proposed 
in [45]; the filter has compact support in a region away from the discontinuities and is optimal 
in the sense that it balances the competing errors caused by localizing it either in physical 
or in momentum space. Its effect in the interacting case is illustrated in figure [181 In this 
logarithmic plot, it extends the range from which one can extract the masses nearly up to 
the midpoint of the lattice. The rather large deviations from the unfiltered data beyond 
that point are irrelevant for our purposes. 

For the two-dimensional models we have extracted the masses from the two point correla- 
tors [30] 

G hos (t) = -l-^^^t + t'.x)^,^)) (97) 

* x V x,x' 



15 In higher dimensions, these Fourier approximation errors away from the discontinuities are negligible in 
comparison to other errors. 
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and 

1 2 

GWm(t) = ^EEE^(* +t >)^'> x ')> ( 98 ) 

* x a=l t> x,x' 

with N t = Ni, N x = AT 2 . Taking spatial averages projects onto = p 2 = 0; this explains 
the observed small discretization errors for the choice r 2 = 4/3 in model (nj, cf. section [4721 
Both correlation functions are proportional to cosh(m(i— N t /2)). The mass can be extracted 
via a linear fit in the logarithmic representation in analogy to the one-dimensional case. 



6 Renormalizibility of the WZ model with SLAC fermions 

Kartens and Smit have demonstrated in [H] that the SLAC derivative induces a nonrenor- 
malizable one-loop diagram in four-dimensional quantum electrodynamics on the latticed 
Therefore, theories involving the SLAC derivative were generally believed to be nonrenor- 
malizable. But it can be shown that the two-dimensional M = 2 Wess-Zumino model is in 
fact renormalizable at least to one-loop order0 

As usual the calculations of lattice perturbation theory are carried out in the thermodynamic 
limit where the number of lattice points tends to infinity and the lattice momentum becomes 
continuous. In momentum space the finite lattice spacing a is translated into a finite cutoff 
A = -. It has to be shown that the diagrams can be renormalized when this cutoff is 
removed. The BPHZ renormalization scheme is used to prove that the renormalized integrals 
tend towards their continuum counterparts, and the counterterms can be identified with 
similar quantities of continuum perturbation theory. The argumentation employed here is 
closely related to the renormalization theorem of Reisz [49]. This theorem does however 
not apply in its original form because the integrands are not smooth functions of the loop 
momentum. 

In the following, we will determine an upper bound for the boson propagator in momentum 
space which will be used later on to argue that parts of the integrals in lattice perturbation 
theory are going to vanish in the continuum limit. For the SLAC derivative, the momentum 

16 In [47], it was suggested that this might be resolved by a different resummation of the perturbation 
series. But in the case of the Schwinger model, this prescription does not have the correct continuum limit, 
as shown in [4*8] . 

17 In fact, the discussion can easily be extended to prove renormalizability also for the Af = 1 model. 
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space representation of the propagators 



1 , -if(k)+m 

and Ju 2 ■ ^ (99) 



P(A;) 2 + m 2 P(£;) 2 + m 2 

for bosons and fermions contains the saw tooth function 

p^k) = - 2/A where (21 - 1)A < < (21 + 1)A. (100) 

The momentum integration is always restricted to the first Brillouin zone, BZ = {(k^) \ \k^\ < 
A}. The internal lines in the one-loop diagrams carry either the internal momentum k or 
a sum k + q of internal and external momenta where q denotes a linear combination of the 
external momenta (using momentum conservation, there are n — 1 such linear combinations 
qj in a diagram with n vertices). Integrations over loop momenta k^ can be split into 
integrations over a square D = {(k^) \ \k^\ < ^} for an arbitrary e < \ and the rest of 
the Brillouin zone, BZ\D. We will argue below that the integral over D converges to the 
continuum value of the integral whereas the integral over BZ\D is shown to vanish as a goes 
to zero. 

Namely, for a given set of external momenta {qj}, one may choose i] = max M j{— |<Zj>|} with 
a small enough such that < rj < e < \. For (k^) e D, we can then read off from 

a\kn ± q^l < a(\k^\ + |<^|) < n(e + rj) for all a < a (101) 

that \k^±q^\ < Ae' with e' :— e + rj < 1, i. e., (k^±q^) is also inside the first Brillouin zone. 
On the other hand, if (k^) e BZ\D, 

n(e -r,)< a(\k^\ - |^|) < a\k» + q„\ < a(\k^\ + |^|) < tt(1 + r,) (102) 

for such lattice spacings a. The latter inequality may be used in order to find an upper 
bound for the propagator, 

< D ,u\^2 < ( 103 ) 



P(k±q) 2 + m 2 P(k±q) 2 
with C = ((e -ri)V2ny 2 . 

It can be easily seen that in the Wess-Zumino model, only two different types of integrals 
contribute at one-loop level. A typical integral of the first type is 

r _fk i 

£ + * ~ J BZ (2vr) 2 (P(kf + m 2 )(P(k + qi ) 2 + m 2 ) . . . (P(k + g n _0 2 + m 2 ) ' ( } 
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d 2 k 

D 



(2tt) 2 (k 2 + m 2 )((k + qi ) 2 + m 2 )... ((k + g n _i) 2 + m 2 ) ' 

d 2 k 1 

BAD (2vr) 2 (k 2 + m 2 )[P{k + qi ) 2 + m 2 )... {P{k + q n ^) 2 + m 2 
d 2 k 1 / 2tt ; 



£ (Co2rl XK^(^0^)^ C ° 2) "" 1Og ( 1 + 



'|*|<^/a (2tt) 2 (A: 2 + m 2 ) V a 2 m 2 , 

Here, we have applied (11021) in order to find an upper bound for the integrand in I n and 
then enlarged the integration domain to a full disk including the first Brillouin zone. Thus, 
I n vanishes in the continuum limit if n > 1. Therefore, the integral I £ tends to the contin- 
uum value of the integral as a goes to zero (and the corresponding continuum integral is 
convergent by power counting), so as long as we are considering diagrams with more than 
one vertex, this type of integrals does not spoil renormalizability. 
Another class of integrals is 



w w f _fk P,{k)P v {k + q 1 )...P e {k + q l ) 

£+ " i BZ (2n) 2 (P(k) 2 + m 2 )(P(k + qi ) 2 + m 2 )...(P(k + < 



IL 



< 



BZ (2tt) 2 (P(k) 2 + m 2 )(P(k + qi ) 2 + m 2 ) . . . (P(k + q n ^) 2 + m 2 ) ' 

d 2 k k^ . . . (k + q~i) e 

(2tt) 2 (k 2 + m 2 )((k + qi ) 2 + m 2 ) . . . ((k + q^) 2 + m 2 ) ' 

d 2 k P M (fc)...P g (fc + g t ) 

BZXD (2tt) 2 (k 2 + m 2 )(P(k + qi ) 2 + m 2 ) . . . (P(k + q n ^) 2 + m 2 ) 
d 2 k \P^k)\...\P s {k + q{)\ 



(105) 



/bz\d (2tt) 2 (k 2 + m 2 )(P(k + Ql ) 2 + m 2 ) . . . (P(k + q n ^) 2 + m 2 ) 

The qi are taken from the qj, so I < n — 1. The same arguments as above show that the 
continuum limit is correct for any n > 2 (again, all continuum integrals are convergent by 
power counting). 

Therefore, renormalizability only remains to be shown for two kinds of integrals. The first 
consists of diagrams with n — 1, e.g., tadpole diagrams. In this case, the loop momentum 
is independent of the (vanishing) exterior momentum so that the argument of P^k) is 
restricted to the first Brillouin zone (where P^{k) = k^). The boundary of the integration 
region behaves as a finite cutoff that is removed in the continuum limit so that the integral 
approaches its continuum counterpart. In the BPHZ renormalization scheme these diagrams 
are just subtracted and do not contribute to the renormalized quantities. 
The second kind of integrals (with n = 2 and I = 1 in (11051) ) requires a more careful investi- 
gation which may be found in appendix [A] This demontrates that lattice discretizations of 
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the two-dimensional Af = 2 Wess-Zumino model based on the SLAC derivative are renor- 
malizable at first order in perturbation theory and yield the correct continuum limit. It 
seems however problematic to use the BPHZ renormalization scheme to renormalize the 
corresponding diagrams in higher-dimensional cases since this would require a differentia- 
tion of the integrands with respect to external momenta. The discontinuity of the saw tooth 
functions in this case would lead to singular terms. 

7 Summary and outlook 

In this paper, we have studied supersymmetric Af = 2 Wess-Zumino models in one and 
two dimensions. The six quantum mechanical models under consideration differ by the 
choice of lattice derivatives and improvement terms. The latter can be used to render the 
theory manifestly supersymmetric on the lattice; in distinction to previous works on this 
subject, our simulations of the broken Ward identities at strong coupling prove that only 
one supersymmetry can be preserved. We have demonstrated to a high numerical precision 
that by far the best results for bosonic and fermionic masses can be obtained from a model 
with Wilson fermions and Stratonovich prescription for the evaluaton of the improvement 
term and from a model based on the SLAC derivative. It is interesting to note that for 
SLAC fermions no improvement term is needed to recover supersymmetry in the continuum 
limit. 

As a key result of this paper for two-dimensional Wess-Zumino models, we propose a non- 
standard Wilson term giving rise to an 0(a 2 ) improved Dirac operator 

o xclt 4 
M = 7 "d M + — 73 A with r 2 = -. (106) 

The masses extracted from this model approach the continuum values much faster than those 
for the model with standard Wilson-Dirac operator (|6T1) . Again, results of a comparably 
good quality can be obtained with nonlocal SLAC fermions. In our case, the common 
reservation that the SLAC derivative leads to non-renormalizable theories (as originally 
shown in [16] for the case of four-dimensional gauge theory) does not hold; we have proven 
that the Wess-Zumino model in two dimensions with this derivative is renormalizable to 
one-loop order. 

Motivated by the fact that the masses for the Af = 2 Wess-Zumino lattice model with 
ultralocal Dirac operator (11061) are quite close to the continuum values already for moderate 
lattices we plan to study the model at strong couplings where we will see deviations from 
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the free theory. We are about to implement the PHMC algorithm |50j as a possibility to 
deal with the small eigenvalues of the fermionic operator. We believe that the Af = 2 
Wess-Zumino model as a simple and well-understood theory without the complications of 
gauge fields has the potential to become a toy-model for developing efficient algorithms for 
systems with dynamical fermions, similar to the ubiquitous Schwinger model which serves 
as toy model for more complex systems with a chiral condensate, instantons, confinement 
and so forth. 

Wess-Zumino models are the flat-space limits of Landau-Ginzburg models. A related project 
might be to study another limit of Landau-Ginzburg theories, namely sigma models in non- 
trivial Kahlerian without a superpotential. Such sigma models admit two supersymmetries; 
typically they have instanton solutions and chiral condensates may be generated. If there 
exist local Nicolai variables which give rise to improved lattice models with one quarter of 
supersymmetry, contact could be made with our investigations of Wess-Zumino models in 
a much broader physical context. Clearly these interesting field theories deserve further 
attention, both from the algebraic and from the numerical side. 

A further obvious problem is to study the nonperturbative sector of the two-dimensional 
Af = 1 Wess-Zumino model. This model shows a richer phase structure than the model 
with two supersymmetries. The sign problem for the Pfaffian seems unavoidable. It is 
interesting to note that in conventions with hermitean gamma matrices, a nonvanishing 
Wilson term for Majorana fermions has to enter the Dirac operator as in f !106p . Nevertheless, 
the discretization errors in this case will be of order 0(a); this can in principle be improved 
to 0(a 2 ) by a slightly less natural substitute for the Wilson term. In general, for the Af = 1 
model, no local Nicolai variables can be constructed which would suggest a supersymmetric 
completion of the naive lattice action. However, one might expect (as in quantum mechanics) 
that such improvement terms for the two-dimensional Af = 1 model with SLAC fermions 
are in fact dispensable. 
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A Renormalization of the fermion loop 



In section M, we have demonstrated which kinds of integrals are potentially dangerous for 
the one-loop renormalizability of the J\f = 2 Wess-Zumino model. The missing integral in 
this proof was given by (11051) with n = 2 and I = 1; it appears for fermion loops with two 
internal lines, 



k 



x / ^ P ^ k ) P ^ ~ g) 

/ BZ {27i) 2 {P 2 {k) + m 2 ){P 2 {k-q)+m 2 )' { ' 



Depending on the superpotential more than two external lines may appear; in this case q 
denotes the sum of all incoming external momenta. In order to facilitate the evaluation of 
this integral, we add (and subtract) a term with a finite continuum limit; apart from that, 
the integral needs to be regulated in this limit. On the lattice with finite lattice spacing, 
BPHZ regularization means to subtract the (as yet finite) value of the integral with vanishing 
exterior momenta. Thus, we consider 

d 2 k P^{k)P^(k-q) f d 2 k m 2 



BZ (2tt) 2 (P 2 (k) + m 2 )(P 2 (k -q) + m 2 ) J BZ (2tt) 2 (P 2 (k) + m 2 )(P 2 (k - q) + m 2 
(value at q — 0) 

d 2 k P^k-q)(k fl -P^k-qj) 



L!Z (2tt) 2 (P 2 (k) + m 2 )(P 2 (k -q) + m 2 ) 

d 2 k q^P^k - q) 

Lsz (2vr) 2 (k 2 + m 2 ){P{k - q) 2 + m 2 ) 



9A V^ r A+q jlc f A dku ^ P^{k-q) 



_ A M J_ A (2tt) 2 (k 2 + m 2 ) (P(k - q) 2 + m 2 ) ' 

In the last step, we have chosen do in such a way that for all a = n/ A < a , shifting k^ £ BZ 
by — g M , one winds up either in the same or in an adjacent Brillouin zone, i. e., 

P^k - q ) = k^-q^ + 2A(6(-A -k^ + q„) - Q(k^ - q„ - A)). (109) 

The first term on the right-hand side of (11081) can be easily seen to converge to the value 
of its continuum counterpart by similar arguments as in (11041) and ( 11051 ). In order to prove 
that the second term does not give rise to any corrections in the continuum limit, we make 
use of (11011 ) and ( 1103H and observe that an upper bound for its modulus is given by 

f A dk ^ \P»(k-q)\ 
-_A M J_ A (27r) 2 (k 2 + m 2 ){P{k — q) 2 + m 2 ) 



2A^ / dk,. 
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< 2CV 

c 
2" 



- A+ql dh 

a 27 
dJ$2 arctan 



A dk 2 1 
_ A 2tt k 2 + m 2 



(qi <-> 92, <-> fe) 



u(k 2 ) - qiA.uj(k 2 



,-1 



A 2 cj(A; 2 



T)u(k 2 



C f A „ 



91 



m 2 + fcf - Aq 1 + A 2 



(gi <-> g 2 ,^i <-> ^2 



+ (<7i <-> 92, fei <-> fc 2 ) 
(110) 



with w(A;) = y/m 2 + A; 2 ; here, we have also used that | arctan(x)| < It is obvious that 
this upper bound converges to zero in the limit where the lattice cutoff is removed. 
This completes the proof that the discretization of the Af = 2 Wess-Zumino model based 
on the SLAC derivative is one-loop renormalizable. 
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